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Strain-engineered A-type antiferromagnetic order in YTiOs: a first-principles 
calculation 
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The epitaxial strain effects on the magnetic ground state of YTiC>3 films grown on LaAlC>3 substrates have 
been studied using the first-principles density-functional theory. With the in-plane compressive strain induced 
by LaAlC>3 (001) substrate, A-type antiferromagnetic order emerges against the original ferromagnetic order. 
This phase transition from ferromagnet to A-type antiferromagnet in YTiC>3 film is robust since the energy 
gain is about 7.64 meV per formula unit despite the Hubbard interaction and modest lattice changes, even 
though the A-type antiferromagnetic order does not exist in any i?Ti03 bulks. 

PACS numbers: 75.25.Dk, 75.30.Kz 
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Transition metal oxides with the perovskite structure 
exhibit a wide variety of electronic phases with plenty 
charge-, magnetic-, and orbital-structures, and show 
many prominent functionalities including colossal mag- 
netoresistance, high-Tc superconductivity, and metal- 
insulator transitions. 1 ' 2 Among these oxides, the i?TiC>3 
family (R is a rare earth cation) , whose 3d ti z bands lay 
near the Fermi level, have not been extensively studied. 
However, i?TiC>3 is not a feature-less family, which also 
owns rich spin/orbital ordered phases.— These phases also 
involve the couplings between charge, orbital, lattice, and 
spin degrees of freedom, which have the potential to be 
used in spintronic or correlated-electron devices. 

i?TiC>3's with trivalent R cations are all protype Mott- 
Hubbard insulators and their common crystal structure 
is a pseudocubic perovskite with an orthorhombic distor- 
tion (the GdFe03-type distortion) . This distortion arises 
from the tilting TiC>6 octahedron around the [110] axis 
and a followed rotation around the [001] axis. The mag- 
nitude of this distortion depends on the ionic radii of R. 
Similar to the _RMn03 case, the lattice structure is more 
distorted with a small R and the Ti-O-Ti bond angle is 
decreased more significantly from 180°. The GdFeC>3- 
type distortion plays a crucial role in controlling the sub- 
tle competitive exchange interactions in these insulating 
titanates. The magnetic ground state of i?Ti03 exhibits 
a transition from ferromagnetic (FM) order to antiferro- 
magnetic (AFM) one with increasing size of R cation4£ 

It is very interesting to compare the phase diagrams of 
i?Ti03 and -RM11O3, both of which show magnetic tran- 
sitions with increasing size of R cations. And the Curie 
(or Neel) temperatures show V-shape behaviors near the 
critical points in both families. However, there are also 
two key differences between _RTiC>3 and i?MnC>3. First, 
the FM-AFM tendency is opposite in these two families. 
With a smaller R, RMn0 3 is more AFM but i?Ti0 3 is 
more FM. Second, the phases revealed in i?MnC>3 are 
more complex than those in i?TiC>3. The AFM phase 
in i?TiC>3 bulks is the simple G-type AFM one while the 
AFM phases (e.g. A-type AFM, spiral-spin order, E-type 



AFM) in _RMn03 are more complex which can be more 
interesting than the simple G-type one£ Thus it is non- 
trivial to ask whether is there any more (hidden) mag- 
netic orders in i?Ti03? In a previous theoretical work 
Ref. 0, total energies of different magnetic structures in- 
cluding A-type AFM, FM and G-type AFM were calcu- 
lated by using an effective spin-pseudospin Hamiltonian, 
which showed that the A-type AFM to FM phase transi- 
tion occurs with increasing GdFe03-type distortion while 
the G-type AFM one has much higher energies. How- 
ever, this result disagrees with the experimental phase 
diagram since the G-type AFM phase is very robust in 
_RTi03 with large R while the A-type AFM phase has 
not been observed in any real i?Ti03 compounds so far. 

In this paper, by using the first-principles calculations, 
we intend to investigate the effects of strain on magnetic 
structures of YTi03 film, focusing on the phase transi- 
tion of the magnetic ground state. Our calculation pre- 
dicts that a robust A-type AFM phase can be stabilized 
by an in-plane compressive strain by using small lattice 
substrates like LaA103. 

YTi03 bulk has a orthorhombic structure (space group 
Pbnm) with lattice constants of a=5.358 A, 6=5.696 A, 
and c=7.637 A. Such a minimum unit cell consists of 4 
formula units. To simulate the effect of in-plane com- 
pressive strain induced by the substrate, the lattice con- 
stants along the a-axis and b-axis are fixed to a=6=5.366 
(3.794 x \/2) A to match the (001) LaA10 3 substrate. 
Here LaA103 is adopted as the substrate to give a weak 
in-plane compressive strain to YTi03 film since the in- 
plane lattice constant of LaA103 is a little smaller (~ 3%) 
than that of YTi03 itself. Such a small difference be- 
tween lattice constants also promises probable epitaxial 
growth of YTiOs thin films on LaA103 substrate. 

Our first-principles calculations were performed us- 
ing density-functional theory (DFT) within the gener- 
alized gradient approximation GGA+U method 8,9 with 
the Perdew-Becke-Erzenhof parametrizatior>i£ as imple- 
mented in the Vienna ab initio simulation package 
(VASP)^i^ The valence states include 4(^ 5s 2 , 3d 2 4p 2 
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TABLE I. The energy difference Ae (per Ti) between mag- 
netic states and the NM state for unstrained bulk YTiOa: 
_E(magnetic)-iJ(NM), in unit of eV and the corresponding lo- 
cal magnetic moments per Ti in unit of ■ 



Magnetic order 


NM FM A-AFM 


C-AFM 


G-AFM 


AE 


-0.533 -0.530 


-0.527 


-0.526 


Magnetic moment 


0.88 0.86 


0.83 


0.82 



and 2s 2 2p 4 for Y, Ti and O, respectively. The lattice op- 
timization and all other static computations have been 
done with the Hubbard U on the (i-electrons of Ti 3+ ion, 
and the Dudarevi^ implementation with U e g — 3.2 eV 
has been used if not noted explicitly^. The atomic posi- 
tions are fully optimized as the Hellman-Feynman forces 
are converged to less than 1.0 meV/A. This optimization 
and the electronic self-consist iterations are performed 
using the plane-wave cutoff of 500 eV and a 9 x 9 x 6 
Monkhorst-Pack /c-point meshi^ centered at T grid in 
combination with the tetrahedron method^ 

First, the ground state of YT1O3 bulk is checked. 
Using the experimental crystal structure, non-magnetic 
(NM) state and four magnetic orders: FM, A-type AFM, 
C-type AFM and G-type AFM, have been calculated to 
compare the energies. Within GGA+U, our calculations 
confirm that the FM order has the lowest energy and 
the calculated local magnetic moment is 0.88/is/per Ti 
in agreement with the experimental magnetic moment 
(0.84/ig)ji£ The detail results of calculated total energy 
are summarized in Table I. According to Table HI other 
magnetic orders' energies (per Ti) are higher than the 
FM one: 3 meV higher for A-type AFM, 6 meV higher 
for C-type AFM and 7 meV higher for G-type AFM. 
It should be noted that the FM ground state is robust 
within a large region of U e s from eV to 5 eV (not shown 
here). Thus, our calculations agree quite well with the 
experimental results and previous DFT studies i 14 i 18 

Subsequently, DFT calculations with the epitaxial 
strain are performed. Epitaxial strain is here realized by 
fixing the in-plane lattice constants to fit the LaA103 sub- 
strate as stated before, while the lattice constant along 
c-axis is varied from 7.0 A to 9.0 A to search the equi- 
librium one under the strain, as shown in Fig. HJa). In 
our calculations, the internal atomic positions are relaxed 
with magnetism under each fixed lattice framework to ob- 
tain optimal crystal structures for calculating accurate 
energies. According to Fig. [Ha), it is obvious that the 
C-type AFM and G-type AFM states are much higher 
in energy than the FM and A-type AFM states. Thus, 
in the following, we will mainly focus on the FM and 
A-type AFM states. The relaxed lattice constant along 
c-axis is 8.25 A for the A-type AFM state and 8.26 A for 
the FM state. These two values are very close, imply- 
ing that the magnetostriction is weak in YTiC>3 , at least 
along the c-axis. And with the optimized c-axis lattice 
constant, the A-type AFM state has a lower energy than 
FM one, e.g. the energy difference between FM and A- 
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FIG. 1. (Color online) (a) Energies for different magnetic 
orders as a function of the c-axis lattice constant, (b) The 
energy difference between the A-type AFM and FM as a func- 
tion of the c-axis lattice constant, (c) The Ti-O-Ti bond angle 
in ab-plane and along the c-axis respectively for the A-type 
AFM state. 



type AFM reaches 7.64 meV per Ti. The A-type AFM 
state appeared in strained YTi03 films is quite nontriv- 
ial since it does not exist in any i?Ti03 bulk, namely 
it is a new phase for i?TiC>3 family. More importantly, 
this strain-induced phase transition from FM to A-type 
AFM is quite promising according to our calculation. As 
shown in Fig. QJb), the energy difference between these 
two orders does not change sign for a large range of o 
axis lattice constant around the optimized one, which 
means this transition is not sensitive to the optimized 
axis lattice constant. Noted that the energy difference is 
relatively significant since the energy difference in bulk 
YTi03 is only 3 meV per Ti. To confirm that this phase 
transition is robust against the change of Hubbard pa- 
rameter, the energy difference between FM and A-type 
AFM states is calculated with different U c s from eV 
to 5 eV stepped by 1 eV, which changes from 16 meV 
to 2 meV (always positive). In other words, this FM to 
A-type AFM transition will not change by varying the 
Hubbard interaction U e g in a large value (from eV to 
5 eV) In short, this strain-induced A-type AFM phase 
will be very promising to be found in real thin films even 
if the experimental lattice constant along c-axis and its 
Hubbard interaction are not exactly the same with those 
in our calculations. 

To understand the underlying physical mechanism, it is 
meaningful to compare the Ti-O-Ti bond angle in YTi03 
with and without the strain, as shown in Table. HU Ac- 
cording to Fig. [Tfc), the bond angle in the a6-plane de- 
creases but the one along c-axis increases with the in- 
creasing c-axis. These results imply that YTi03 is com- 
pressed and thus more distorted in the a6-plane but elon- 
gates and thus is less distorted along the c-axis. 

As stated before, it is well known that in i?Ti03 
compounds, small Ti-O-Ti bond angles with more dis- 
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TABLE II. Bond angles in the afe-plane and along c-axis of 
YTiQ 3 film on LaA10 3 substrate and bulk YTi0 3 . 



Ti-O-Ti 


YTi0 3 film 


YTi0 3 film 


bulk YTi0 3 


bond angle 


(A-AFM) 


(FM) 


(FM) 


a6-planc 


139.6° 


139.5° 


144.3° 


c-axis 


143.7° 


143.5° 


141.9° 
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FIG. 2. (Color online) (a)Total DOS of YTiOs film. The 
Fermi energy is positioned at zero, (b) The energy gap as a 
function of U e ff. The critical U g value for zero gap is esti- 
mated as 0.6 eV by extrapolation. 



torted lattice structure exhibit FM order while those 
with larger Ti-O-Ti bond angle tend to be AFM. It 
should be noted that this tendency is in opposite to the 
tendency in manganites. In manganites, the straighter 
(closer to 180°) bond angle should lead to ferromagnetic 
exchange. Therefore, the strain driven phase transition 
here is not exactly the same with the corresponding one 
in strained manganitesJ^ So here in strained case, the 
spin order in a6-plane has a tendency towards FM corre- 
lation due to decreasing in-plane Ti-O-Ti bond angle, 
while the increasing bond angle along c-axis tends to 
drive spins arranged antiparallel. In this sense, the emer- 
gence of A-type AFM in YTi0 3 on LaA10 3 substrate 
can be qualitatively understood as the ferromagnetism 
with decreasing bond angle in afe-plane and antifcrro- 
magnetism with increasing bond angle along c-axis. Of 
course, a comprehensive understanding of this magnetic 
transition needs more careful studies from microscopic 
theory. Furthermore, the bond angles are very close be- 
tween FM and A-AFM states in strained films, implying 
that the exchange-striction effect is weak in this materi- 
als, which is different from the strong exchange-striction 
in manganites^ or ion-selenidesi^i 

In many correlated electron materials, accompanying 
magnetic transitions, conductance often changes drasti- 
cally, e.g. metal-insulator or insulator-superconductor 
transitions. Therefore, it is necessary to check the den- 
sity of states (DOS) of YTi03 under strain, as shown in 
Fig. Ha). 

First of all, YTiOs bulk is an insulator in our DFT 
calculation (not shown here), in agreement with the ex- 



perimental result 4 This insulating behavior is due to 
the Coulomb repulsion between 3d electrons, implying 
a Mott insulator. This scenario can be easily demon- 
strated because the pure GGA calculation without the 
Hubbard U gives a metallic DOS. The band gap is 1.504 
eV for bulk YTiOs, a little overestimated compared with 
the experimental data 1.2 eV£2- Our calculation found a 
gap as 1.564 eV for YTi03 under strain, which is only a 
little larger than the bulk value. This gap exists within 
a wide range of U e g from 1 eV to 7 eV, though the detail 
value depends on U e s, as shown in Fig. &[b). Therefore, 
it is safe to say that the strained YTiOs remains an insu- 
lator and its conductance is not obviously changed since 
the gaps are almost identical. 

In summary we have studied the effects of epitaxial 
strain on the magnetic ground states in YTiOs films. 
Our results predicted a new magnetic ground phase A- 
type antiferromagnet which had not been realized in any 
i?Ti0 3 bulk compounds. This robust A-type AFM phase 
is stabilized by an in-plane compressive strain induced 
by LaA103 substrate. Its origin is understood as the fer- 
romagnetism with decreasing bond angle in the a&-plane 
and antiferromagnetism with increasing bond angle along 
the c-axis. Furthermore, the density of states calcula- 
tion confirmed that the insulating behavior and the en- 
ergy gap would not be significantly affected by this strain 
driven magnetic transition. 
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